SELFDUAL SUBSTITUTIONS IN DIMENSION ONE 



VALERIE BERTHE, DIRK FRETTLOH, AND VICTOR SIRVENT 



Abstract. There are several notions of the 'dual' of a word/tile substitution. We show that 
the most common ones are equivalent for substitutions in dimension one, where we restrict 
ourselves to the case of two letters/tiles. Furthermore, we obtain necessary and sufficient 
arithmetic conditions for substitutions being selfdual in this case. Since many connections 
between the different notions of word/tile substitution are discussed, this paper may also 
serve as a survey paper on this topic. 



1. Introduction 

Substitutions are simple but powerful tools to generate a large number of nonperiodic struc- 
tures with a high degree of order. Examples include infinite words (e.g., the Thue Morse 
sequence, see [Ml EH]), infinite tilings (e.g., Penrose tilings, see [381 03]) an d discrete point 
sets (e.g., models of atomic positions in quasicrystals, see [31 [31]). Here we consider several 
instances of the concept of substitutions: 

(a) word substitutions 

(b) nonnegative endomorphisms of the free group Fi = (a, b) 

(c) tile-substitutions 

(d) dual maps of substitutions 

Each of the concepts above gives rise to the concept of a dual substitution. Our first goal is 
to show the full equivalence of the distinct notions of dual substitution with respect to the 
concepts above for the case of two letters, or of two tiles in R 1 . Thus we will exclusively 
study substitutions on two letters, and for the sake of clarity, we will define every term for 
this special case only. From here on, we will consider only this case, whether or not this is 
stated explicitly. 

Let us mention that there is a wealth of results for word substitutions on two letters, thus for 
tile-substitutions in R 1 with two tiles. For a start see |36 } 1291 [35] and references therein. The 
following theorem lists some interesting results which emerged in the work of many authors 
during the last decades. We recall that a word substitution is a morphism of the free monoid. 
A word substitution is said to be primitive if there exists a power n such that the image of 
any letter by a n contains all the letters of the alphabet. 

Theorem 1.1. Let a be a primitive word substitution on two letters. Then the following are 
equivalent: 
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(1) Each bi-infinite word u generated by a is Sturmian, i.e., u contains exactly n + 1 
different words of length n for all n G N. 

(2) The endomorphism a : F2 — > is invertible, i.e., a G Aut(F2). 

(3) Each tiling generated by a is a cut and project tiling whose window is an interval. 

(4) There is k > 1 such that the substitution a k is conjugate to L ai EL a2 EL a3 ■ ■ ■ EL am , 
where G N U {0} , L : o — >■ a, 6 — > ab, E : a — > b, b — > a. 

The equivalence of (2) and (4) follows from [45J. The equivalence of (2) and (3) appears in 
|26| . and relies on earlier results, see references in [26, [H]. The equivalence of (1) and (3) is 
due to [321 E]. For the equivalence of (1) and (4), see [30]. For more details, see [29| [35]. 

This theorem links automorphisms of the free group F% with bi-infinite words over two letters, 
and with tilings of the line by intervals of two distinct lengths. Our motivation was to work 
out this correspondence. 

The paper is organised as follows: Section [2] provides the necessary definitions and facts about 
word substitutions, and about (inverses of) substitutions as automorphisms of F2. Section [3] 
contains the relevant definitions and facts about (duals of) substitution tilings and cut and 
project tilings. The very close relation of the natural decomposition method of [32] to duals of 
cut and project tilings via so-called Rauzy fractals is given in Section I3~3l Section |4]introduces 
generalised substitutions and the essential definitions and facts for strand spaces and dual 
maps of substitutions. Any reader familiar with any of this may skip the corresponding 
section(s). Section [5] explains how the distinct concepts of dual substitutions are related. The 
first main result, Theorem 15.101 shows that they all are equivalent, in the sense, that they 
generate equivalent hulls. (For the definition of the hull of a word substitution see Equation 
(|2|) below. For a precise definition of equivalence of hulls, see Section [5j) In Section [6] we give 
several necessary and sufficient conditions for a substitution to be selfdual, namely, Theorems 
16.41 and 16.71 These results yield a complete classification of selfdual word substitutions on two 
letters. 

2. Word substitutions 

2.1. Basic definitions and notation. Let A = {a, b} be a finite alphabet, let A* be the 
set of all finite words over A, and let A 1, = {(tij)iez I u i £ A} be the set of all bi-infinite 
words over A. A substitution a is a map from A to A* \ {e}; it is called non-erasing if it 
maps A to A* \ {e}, where e denotes the empty word. We will exclusively consider here 
non-erasing substitutions. By concatenation, a substitution extends to a map from A* to ^4*, 
and to a map from A z to A 1 ". Thus a substitution can be iterated. For instance, consider the 
Fibonacci substitution 

(1) a-.A^-A*, a (a) = ab, a(b) = a. 

We abbreviate this long notation by a : a — >• ab, b — >■ a. Then, a{a) = ab, <r 2 (a) = a (a (a)) = 
a(ab) = a(a)o~(b) = aba, <r 3 (a) = abaab, and so on. In order to rule out certain non-interesting 
cases, for instance E: a 1— >■ b,b 1— >■ a, the following definition is useful. 

Definition 2.1. Let a be a substitution defined on A. The substitution matrix M a associated 
with a is defined as its Abelianisation, i.e., 



M <y = {\°-(j)\i)i<i,j<2, 
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where \w\i stands for the number of occurrences of the letter i in w. 

A substitution a on A is called primitive, if its substitution matrix is primitive, that is, if 
(M a ) n > for some n G N. 

A substitution a is called unimodular, if its substitution matrix M a has determinant 1 or — 1. 

In the sequel, we will consider unimodular primitive substitutions only. The requirement of 
primitivity is a common one, essentially this rules out some pathological cases. The require- 
ment of unimodularity is a restriction which we use because we will focus on the case where 
the substitution is invertible (see Section [2.21 for a definition). 

We will consider the set of all bi-infinite words which are 'legal' with respect to a substitution 
a. One reason for this choice of working with bi-infinite words (instead of infinite words) 
comes from the connections with tile-substitutions that we stress in Section [3l Thus we 
define the hull X a of a as follows. 

(2) X a = {u G A z | each subword of u is a subword of cr n (a) or a n (b) for some n}. 
Subwords of a word are commonly also called 'factors' of that word. 

Proposition 2.2. |36|, I35| For any primitive substitution a and for alln > 1 holds: X a ™ = X a . 

There is a further notion of a hull, namely, the hull of a bi-infinite word. In our framework 
(substitutions are assumed to be primitive) the two notions coincide. The hull X u of a bi- 
infinite word u is defined as the closure of the orbit {S k u \ k G Z} in the obvious topology, 
see [361 E5J- Here S denotes the shift operator Su = S{ui)i & % = (wj + i)j 6 z, which yields 
S~ x u = («i_i)i 6Z . 

Proposition 2.3. Let a be a primitive substitution. Then, for each u G X a , X a = X u . In 

particular, X a is determined by each u G X a uniquely. 

Let a be a primitive substitution. By the Perron-Frobenius theorem, its substitution matrix 
M a has a dominant real eigenvalue A > 1, which has a positive eigenvector (i.e., an eigenvector 
having all its components > 0). We call A the inflation factor of a (see Section 13.11 for a 
justification of this term). If (1 — a, a) is the normed eigenvector of A (0 < a < 1), that is, 
the vector whose sum of coordinates equals 1, then a is called the frequency of the substitution. 
Indeed, it is not hard to see that every bi-infinite word in X a has well defined letter frequencies 
1 — a and a, see for instance [36]. Left eigenvectors associated with the inflation factor A will 
also play an important role here. We will work with the positive left eigenvector v\ = 
with first coordinate 1. 

In all that follows, cr{i)k stands for the fc-th letter of cr(i), \o~(i)\ is the length of a(i) and 
a(i)[k — 1] is the prefix of length k — 1 of o~(i) (with a(i)o being the empty word), i.e., 

o-(i) = a {i)i ■ ■ ■ 0-(«)[<T(i)| = cr(i)[k - l]a(i) k a(i) k+ i ■ ■ ■ ^{i)\ a {i)\- 

We denote by A the Abelianisation map from A* to Z 2 : if w is a word in A*, then A(w) is 
the vector that counts the number of occurrences of each letter in w, i.e., A: A* — > Z 2 , w i— )■ 
(|u;| a , \w\b)- There is an obvious commutative diagram, where M a stands for the substitution 
matrix of cr: 
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A* 

I? 



a 



A* 
| A 

I?. 




2.2. Substitutions as endomorphisms of F2. Naturally, any word substitution on A 
gives rise to an endomorphism a of F2, the free group on two letters. (Note that not 
every endomorphism of F2 gives rise to a proper word substitution, consider for instance 
a : a — > ab^ 1 , b — > b .) If this endomorphism a happens to be an automorphism, that 
is, a £ Aut(F2), then a is said to be invertible. It is well-known that the set of invertible 
word substitutions on a two-letter alphabet is a finitely generated monoid, with one set of 
generators being 

(3) E : a i-> b, b i-> a, L : a h-> a, b h-> ab, L: a i-> a, b 1— )• 6a. 

For references, see [45] and Chap. 2 in |29j . 

In the sequel, we want to consider whether two substitutions generate the same sequences, 
that is, the same hull. Recall that an automorphism 7 is called inner automorphism, if there 
exists w £ F2 such that 7(2;) = wxw~ l for every x £ F?. We let j w denote this inner 
automorphism. 

Definition 2.4. We say that two given substitutions a and g are conjugate, if a = 7^ o g for 
some w £ F2. In this case, we will write a ~ g. 

Remark 2.5. We use here the term 'conjugate' but a more precise terminology would be the 
following: a and g belong to the same outer class, or else, g is obtained from a by action of 
an inner automorphism (also called conjugation). For convenience, we take the freedom to 
use the short version. Note that if a and g are conjugate, then w is a suffix or a prefix of g(x) 
for x = a,b, and o~(x) has the same length as the word g{x) for every x £ A. 

Recall the following theorem of Nielsen |34| : given two automorphisms a and g of the free 
group F2, they have the same substitution matrix if and only if they are conjugate. (Here, the 
substitution matrix relies on the Abelianisation map from F2 to Z 2 that counts the number 
of occurrences of each 'positive' letter minus the number of occurrences of each 'negative' 
letter.) We thus conclude in terms of substitutions (see [3D] for a combinatorial proof): 

Theorem 2.6. Two invertible substitutions a and g are conjugate if and only if they have 
the same substitution matrix. 

2.3. A rigidity result. The following theorem is a classical rigidity result for two-letter 
primitive substitutions by P. Seebold [3D], see also [25] . Here we provide two alternative 
proofs. Let us note that the apparent simplicity of the first proof relies on the use of Theorem 
12.61 and the defect theorem (Theorem 1.2.5 in [28J), whereas the original proof of [40J was 
self-contained. 

Theorem 2.7. Let a, g be primitive substitutions on the two-letter alphabet A. If o~ k ~ g m 
for some k, m, then X a = X Q . 

Furthermore, if a and g are invertible, then X a = X Q if and only if a k ~ g m for some k,m. 
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In plain words, this theorem states that two invertible substitutions are conjugate (up to 
powers) if and only if their hulls are equal. In even different words (compare Proposition 
\2.3\i : if u is a bi-infinite word obtained by some invertible primitive substitution on two 
letters, where the substitution is unknown, then u determines the substitution uniquely, up 
to conjugation and up to powers of the inflation factor. An immediate consequence is the 
following result. 

Corollary 2.8. Let a, g be primitive invertible substitutions on the alphabet A = {a, b} with 
the same inflation factor. Then, a ~ g if and only if X a = X g , which is also equivalent to a 
and g having the same frequency a. □ 

Proof. Let us prove Theorem 12.71 Because of Proposition 12.21 we can restrict ourselves to the 
case k = m = 1. First we prove that if a ~ g, then X a = X g . Thus let a ~ g. This means 
there is w G F% such that a(x) = wg(x)w~ 1 for x = a,b. Consequently, for x G A holds: 
wg(x) = a(x)w, and in general for u G A*, wg{u) = a(u)w, which yields for x G A, 

(4) wg(w) ■ ■ ■ g k ~ 1 (w)g k (x) = a k (x)wg{w) ■ ■ ■ g k ~ 1 {w). 
Setting k;W := wg(w) ■ ■ ■ g k ~ 1 (w), we can write shortly 

w ik) g k {x) = a k (x)w {k \ 

Let x = a. By the defect theorem (Theorem 1.2.5 in [28]) this commutation relation implies 

(5) g k (a) = u k v k , a k (a) = v k u k 

for some words u k and v k , with u k being nonempty. Now, let u G X a , and let v be some 
subword of u. Then, by primitivity, v is a subword of a k (a) for some k. Since m&x(\u k \ , | v k \) — > 
oo for k — > oo {a is primitive), v is also contained in u k or v k for k large enough. Thus it is 
also contained in g k (a) for some k, thus v is a word in X g . The same argument holds vice 
versa, thus X a = X g . 

For the other direction, we first note that if A is an irrational eigenvalue of a 2 x 2 integer 
matrix M then it is an algebraic integer, and the algebraic conjugate A' of A is the second 
eigenvalue of M. If v = (l,v\) is an eigenvector of M corresponding to A, then v' = (l,f\) 
(again, v' x denotes the algebraic conjugate of v\) is an eigenvector corresponding to A'. If the 
matrix is furthermore primitive, both vectors are distinct by the Perron-Frobenius theorem, 
and thus linearly independent. 

Let a, g be two primitive and invertible substitutions such that X a = X Q . Let u G X a . Let us 
recall that the the vector of frequencies (1 — a, a) of letters in u is an eigenvector of the (up 
to here unknown) substitution matrices M a and M g . Hence the vectors (1,^) = (1, jr^) and 
(1,£ ; ) are eigenvectors for both matrices. We now consider the eigenvalues associated with 
the previous eigenvectors. Since a is invertible, M a is unimodular. Its (up to here unknown) 
inflation factor A is therefore a unit in the underlying ring of integers of the form Z[vfc] for 
some k > 1. It is well known that the unit group of Z[vfc] is generated by a fundamental unit 
z. (This is a consequence of the fact that there is a fundamental unit for the solution of the 
corresponding Pell's equation, or a consequence of Dirichlet's unit theorem, see for instance 
|27j or [33]). Thus A is a power of the generating element z. Let A = z n , where n is arbitrary 
but fixed. The same holds for the inflation factor of M e which also belongs to and 
which is thus of the form z m . By algebraic conjugation we obtain the second eigenvalue, and 
the second eigenvector of M g . 
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w (k) p k (a) w (k) p k (a) w (k) p k (a) 
1 1 " 1 " h \~ 1 



a k (a) w (k) o k (a) w (k) a k (a) w (k) 

Figure 1 . Three cases in the proof of Theorem 12.71 

Since the eigenvectors and (1,/) are linearly independent eigenvectors, on which the 

substitution matrices of a n and g m act in the same way, a n and g m have thus the same 
substitution matrix. By the fact that all invertible substitutions with the same substitution 
matrix are conjugate (see Theorem 12 .6p . the claim of the theorem follows. □ 

Remark 2.9. One can prove also directly from (JJD that g k (x) and a k (x) share a common 
subword, without using the defect theorem. This can be done as follows: We prove Equation 
((5|) by the following argument: 

We have to distinguish three cases (compare Figure [1]) . 

Case 1: If \w^\ = \g k {a)\ = \a k (a)\, then = g k (a) = a k (a) (Figured! left). 

Case 2: If \w^\ < \g k (a)\ = \cr k (a)\ (Figure [[J centre), then is a prefix of g k (a) and 
a suffix of a k (a). Moreover, the remaining suffix of g k (a) overlaps the remaining prefix of 
a k (a). Thus we obtain again Equation ([5|): 

(6) g k (a) = u k v k , a k (a) = v k u k 

for some nonempty words u k {= w^) and v k . 

Case 3:: If \w^\ > \g k (a)\ = \cr k (a)\ (Figured! right), then g k (a) is a suffix of w^. Either 
this suffix of overlaps already with cr k (a) (as in the figure), then (|5|) holds for some 
nonempty words u k ,v k . Or (if |u/ fe )(a)| > 2\g k (a)\) g k (a)g k (a) is a suffix of w^ k \ and we 
continue with the shorter words vs a k (a)w^ without the suffix g k (a). After finitely 
many steps, we are in one of the first two cases. 

In either case, ([5]) holds for some nonempty word u k . (In Case 1, just let v k be the empty 
word.) 

Remark 2.10. Let us note that the assumption that a is invertible is crucial in Theorem 12. 71 
as shown by the following example (see [25]). Consider on the alphabet {a, b} the following 
two primitive substitutions: 

a : a i-> ab, b i-> baabbaabbaabba, g : a h-> abbaab, b h-» baabbaabba. 

One has a(ab) = g(ab) and a(ba) = g(ba). We deduce that a and g have the same fixed 
point beginning by a, and thus X a = X g . Nevertheless, a and g are neither conjugate, nor 
conjugate up to a power of a common substitution (their substitution matrices are neither 
conjugate in GL(2,Z), nor conjugate up to a power to a common matrix). 



3. Tile-substitutions 
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3.1. Substitution tilings. In contrast to word substitutions, which act on symbolic objects 
like words, tile-substitutions act on geometric objects, like tiles or tilings. A tiling of M. d is a 
collection of compact sets which cover topologically M d in a non-overlapping way, that is, the 
interiors of the tiles are pairwise disjoint. In R 1 there is a natural correspondence between 
bi-infinite words and tilings when the tiles are intervals: just assign to each letter an interval 
of specified length. 

In general, a tile-substitution in M d is given by a set of prototiles T±, . . . ,T m C M d , an ex- 
panding map and a rule how to dissect each expanded prototile into isometric copies of some 
prototiles T^. Here, we restrict ourselves to two prototiles in dimension one, and, moreover, 
our prototiles are always intervals, with the expanding map being given by an inflation factor 
A > 1. (For the discussion of analogues of some results of the present paper for the case of 
more general tilings, see [38l H31 [18] ) . The precise definition of a tile-substitution in M, where 
the prototiles are intervals, goes as follows. 

Definition 3.1. A (self-similar) tile-substitution in R is defined via a set of intervals T\, . . . T m 
— the prototiles — and a map s. Let 

m 

(7) XT j = (jT i + V ij (l<i<m), 

i=i 

where the union is not overlapping (i.e., the interiors of the tiles in the union are pairwise 
disjoint), and each T>- is a finite (possibly empty) subset ofM. d , called digit set. Then 

afy) := {T i +V lj \i = l...m} 

is called a tile-substitution. It is called primitive if the substitution matrix M s := (|2^jj|)i<ij<2 
is primitive, where \Dij\ stands for the cardinality of the set T>ij. 

By s(Tj + x) := s(T) + Xx and s({T, T'}) := {s(T), s(T')}, s extends in a natural way to all 
finite or infinite sets of copies of the prototiles. 

In analogy to word substitutions we want to deal with the space X s of all substitution tilings 
arising from a given tile-substitution. Note the correspondence with the definition of the 
hull of a word substitution, see Equation ([2]). The main difference between both types of 
associated dynamical systems is that X s is endowed with a Z-action by the shift, whereas X s 
is endowed with an K- action defined by the action of translations. 

Definition 3.2. Let s be a primitive tile- substitution with prototiles T\,T<i- The tiling space 
X s is the set of all tilings T, such that each finite set of tiles ofT is contained in some translate 
of s n (Ti) or s n {T2). Any element o/X s is called a substitution tiling generated by s. 

Remark 3.3. Any primitive self-similar tile-substitution is uniquely defined by its digit set 
matrix T>. This holds because one can derive the inflation factor A and the prototiles Tj from 
the digit set matrix T> := (T>ij}... This is not only true for two tiles in one dimension, but for 

any self-similar tile-substitution in W 1 . For details, see [18]. Here we just mention two facts: 
the inflation factor A is the largest eigenvalue of the primitive substitution matrix M s . And 
the prototiles are the unique compact solution of the multi component IFS (iterated function 
system) in the sense of [30] (also called graph-directed IFS), which is obtained by dividing ([7]) 
by A. In particular the vector of lengths of the tiles is a left eigenvector of the substitution 
matrix M a . 
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Remark 3.4. A one-dimensional tile-substitution (where the tiles are intervals) yields a 
unique word substitution: just replace the tiles by symbols. Conversely, one can realize any 
primitive word substitution as a tile-substitution by taking as lengths l a ,h for the prototiles 
the coordinates of a positive left eigenvector associated with its inflation factor A: the intervals 
Tj are chosen so that they line up with the action of the word substitution. We chose here to 
normalise the eigenvector v\ = (1,^a) = (la,h) by taking its first coordinate equal to 1. For 

j = a,b, let Tj = [0,lj]. For j = a,b, if a(j) = cr(j)i • • • <r(j)\a(j)\ then Xlj = EL=i l cr(j) k , i-e., 
the tile Tj is inflated by the factor A, so it can be subdivided into translates of the prototiles 
according to the substitution rule: 

XTj = [0, Xlj] ^ {T a{j)l , T a{jh + l a(j)l T a(J)wm + l a(j)l +■■■ + la{j) wm _-}- 

This can be written as 

**i= U T j + S([a(j)} k ^), 
hi ■ (j,k)eFi 

where 

F = {(j,k)\j eA,l<k<\a(j)\} 

and 

F t = {(j,k) £ F\a(j) k =i}, 
with the valuation map 5 : A* — > M + being defined for any w = w\ ■ ■ ■ w m £ A* as 

6(wi . ..w m ) = l Wl H 1- l Wm = \ w\ a + \w\ b £ x = (A(w),v x ), 

by recalling that A stands for the Abelianisation map. We then set for all i,j 

Vij = {6([v(j)] k -l) I (j,k)€Fi}. 

We illustrate this by the following example. 

Example 3.5. Consider the square of the Fibonacci substitution from Equation (UJ, namely, 
g = a 2 : a — > aba, b — > ab. We work here with the square of the Fibonacci substitution since 
the determinant of its substitution matrix equals 1. We will realize it as a tile-substitution 
as follows (see Figure EJ): Let T x = [0, 1], T 2 = [0, 1/r], where r = ^±1. Note that A = r 2 is 
the dominant eigenvalue for M (T 2, with (1, 1/r) being an associated left eigenvector. Then 

r 2 T x = [0, 2 + 1/r] = T x U (T 2 + 1) U (T x + 1 + 1/r); r 2 T 2 = [0, 1 + 1/r] = T x U (T 2 + 1), 

where the unions are disjoint in measure. (Note that r = 1 + 1/r.) Hence the last equation 
yields a tile-substitution s: 

s(Tx) = {T U T 2 + l,Tx +r}, s(T 2 ) = {Ti,T 2 + 1}. 

For an illustration of this substitution, see Figure [2j This substitution s can be encoded in the 
digit sets D^i = {0,r},2? 2) i = {1},2?i j2 = {0},P 2j2 = {1}. This can be written conveniently 
as a digit set matrix: 

{0,r} {0y 
{1} {1}, 

By comparison with Definition 12.11 we note that we can derive the substitution matrix from 
the digit set matrix simply as follows: M a 2 = D = (|2\,-|)i<i.j< 2 . In this case we get the 
matrix ( \ \ ) . Its dominant eigenvalue is the inflation factor A = r 2 . 



(8) V 
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Figure 2. The Fibonacci tile-substitution s (box top left), some iterates of 
s on T%, and the generation of a Fibonacci tiling as a cut and project tiling 
(right). The interval W defines a horizontal strip, all lattice points within this 
strip are projected down to the line. 



3.2. Cut and project tilings. In order to define a notion of duality for tile-substitutions, 
we now work in the framework of cut and project sets. 

Certain substitution tilings can be obtained by a cut and project method. There is a large 
number of results about such cut and project sets, or model sets, see [3T] and references 
therein. In our setting, this is pretty simple to explain, compare Figure [2j Let G = H = R, 
let tt\ : G x H — > G, 7T2 : G x H — > H be the canonical projections, and let A be a lattice 
in G x H = R 2 , such that tt\ : A — > G is one-to-one, and ^(A) is dense in H. Then, choose 
some compact set W C H with W being the closure of its interior, and let 

V = {ki(x) I x € A, 7r 2 (x) G W}. 

Then V is a cut and project set (or model set). Since V is a discrete point set in R = G 
(W is bounded), it induces a partition of R into intervals. Regarding these (closed) intervals 
as tiles yields a tiling of R. Such a tiling is called cut and project tiling. Furthermore, a 
tile-substitution s is said to yield the cut and project set V if V is the set of left endpoints of 
an element of X s . According to Theorem II. 1\ a primitive word substitution a on two letters 
yields cut and project tilings whose window W is an interval if and only if it is invertible |26j . 

Given a primitive tile-substitution s which is known to yield a cut and project set V, one can 
construct A and W out of s in a standard way. In general, A and W are not unique. The 
following construction provides a choice of A and W when s comes from a primitive two-letter 
substitution, according to Remark [331 It has the advantage that everything can be expressed 
in some algebraic number field, which allows the use of algebraic tools. See also [20] for more 
details on cut and project schemes in the case of G = H = R. 

We thus start with a primitive and unimodular word substitution a on a two-letter alphabet 
such that its associated tile-substitution s yields the cut and project set V. Since the substi- 
tution M a is an integer matrix, its eigenvalues A and A' are two conjugate quadratic irrational 
numbers. (The case where A is an integer requires that the internal space H is non-Euclidean 
[HE]. Since our substitutions will always be unimodular in the sequel, this case cannot occur 
here.) Let v\ = be the left eigenvector associated with the dominant eigenvalue A. This 

eigenvector yields the 'natural' lengths of the prototiles. Thus let T\ = [0,1], T 2 = [0,£a]- 
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Note that 1 and £ are rationally independent. Now, let 

(9) A = (v,w)z = {«(}) \ a,PeZ}, 

where £' x denotes the algebraic conjugate of £\. The projections iri and 1T2 correspond to the 
canonical projections. Now, consider the set V of endpoints of the intervals in a tiling in X s . 
Without loss of generality, let one endpoint be 0. Then all other endpoints are of the form 
q + j3£ x G %[£\\ C Q(A). Any point a + (3£\ G V has a unique preimage in A (since 1, £\ are 

rationally independent) , namely, a ( } ) + ^ ) . Thus, each point a + j3£\ in V has a unique 

counterpart in the internal space H, namely tt2 o 71 J" (a + /3l\) = a + (3£' x . The map 712 o ir^ 1 
is called star map and will be abbreviated by *. 

Since 712(A) is dense in H and is the closure of its interior, the window W is obtained as the 
closure of 7T2 o tt^ (V). Note that the fact that V is a cut and project tiling guarantees that 
W is a bounded set, since being compact. With the help of the star map we can write shortly 
W = cl(V*). In our context, the star map has a very simple interpretation: by construction 
of the lattice A, the star map is just mapping an element of Q(A) to its algebraic conjugate 

(10) (a + l3£ x )* = a + f3£' x . 

Nevertheless, in general it is more complicated, and one should keep in mind that the star 
map maps G to H. In general, G and H can be very different from each other, for instance 
of different dimension. 

Now recall that any self-similar tile-substitution is uniquely defined by its digit set matrix T> 
(Remark I3.3|) . This allows us to define the star-dual of a tile-substitution by applying the 
star map to the transpose T> T of the matrix P, as performed in [44U19t fT8] where this 'Galois' 
duality for tile-substitutions is developed. 

Definition 3.6. Let s be a primitive invertible self-similar tile- substitution yielding cut and 
project tilings, with digit set matrix T> . Then the star-dual substitution s* of s is the unique 
tile- substitution defined by (V T y , with the star map being defined in (Q2p. 

Here X* means the application of the star map to each element of some set X C Q(A) 
separately. This definition together with Definition 13.21 defines the star- dual tiling space X s *. 

Example 3.7. The star-dual of the tile-substitution associated with the squared Fibonacci 
substitution in Example 13.51 is easily obtained by applying the star map to the transpose of 
the digit set matrix in (|8|). We obtain 

{0} {1}J ' 

3.3. Natural decomposition method and Rauzy fractals. We now explain how to as- 
sociate a candidate window W with any primitive two-letter substitution a in such a way 
that the corresponding tile-substitution s yields a cut and project set. This candidate is the 
so-called Rauzy fractal associated with the substitution a that is introduced below. 

In this section we follow |42j . originally defined on a <i-letter alphabet. We restrict here to the 
case d = 2. As above, let a be a primitive unimodular substitution on the two-letter alphabet 
A = {a,b}. Let A > 1 be its inflation factor and A' the other eigenvalue of the substitution 
matrix M a . Let u = (u n )n£N be an infinite word such that a(u) = u. (Note that we use 
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one-sided infinite words here, in accordance with [32]. The extension to bi-infinite words is 
straightforward.) Since a is primitive, it suffices to replace a by a suitable power of a for 
such a fixed point word to exist. 

The set X a is mapped into R via a valuation (compare with the valuation map 5 introduced 
in Remark I3.4p . which is a map A : A* — > R, satisfying A(vw) = A(v) + A(w) and A(o~(w)) = 
X'A(w) for all v,w £ A*. Hence, A satisfies for all w G A* 

A(w) = \w\ a + \w\bh' = (A( w )i v \')> 

where vy = (l,£y) is the normalised left eigenvector associated with the eigenvalue A' of the 
substitution matrix M a . According to |37j . the set 

9\ := {A(uqUi . . . u m ) | m > 0} 

is called the Rauzy fractal associated with the one-sided fixed point u of the substitution a. 
Since a is unimodular, |A'| < 1 and one deduces that 9\ is a compact set. For more about 
Rauzy fractals, see [35], [8] or else [10] and the references therein. 

Let 

%, := {A(u ui . . . u m ) | u m +i =i, m > 0} 
where i E A. Clearly D\ = £H a U 9V We shall call (9^ a ,!tRfc) the natural decomposition of 
the Rauzy fractal 9^. The natural decomposition of 9\ is the attractor of a graph directed 
IFS (cf. [3U] and also [22]) in the following way, with the notation of Remark 13.41 (D\ a ,^-b) 
satisfies 

(11) «i= |J (\% + A{[a(j)]k-i), for i£A, 

with Fi being defined in Remark 13.41 To prove it, we use the fact that o~(u) = u. 

Note that the sets SHj are the closure of their interior and their boundary has zero measure, 
as proved in [42] in general case of a d-letter substitution. For more properties, see e.g. 
\35\ l8| [TO]. Furthermore, the sets 9^ are not necessarily intervals. In fact it can be proved 
that they are intervals if and only if the two-letter primitive substitution a is invertible (see 
also Theorem 1 1.1 j) . For more details, see e.g. |26 } 113 } [Tl] for proofs of this folklore result. 

A moment of thought yields that A maps the m-th letter of u, which corresponds to the right 
endpoint |u| a + of the m-th tile in the tiling in G (corresponding to the word u) to 

\u\a + M&^a i n H- I n other words, 9\ is a right candidate for the window W for the tiling in 
G, with V = S({uq . . . u m | m > 0}) (see Remark 13.41 for the definition of the map S): 

%\ = 7T2 O 7Pf {6(UQ ■ ■ ■ Um) \ TU > 0}, 

with 7Ti and 7T2 being defined in Section [3721 In particular 7T207r^ 1 (5(uo • • • u m) = A(uq . . . u m ), 
for any m > 0. It remains to check that {n\{x) \ x € A, ^(x) € = V, i.e., that we do not 
have V strictly included in {tti(x) \ x € A, ^(x) G This comes from the following result. 

Theorem 3.8. [5j [23| [21] Let a be a primitive unimodular substitution on two letters. Then 
a yields the cut and project set V = {8(uq . . . u m ) \ m > 0} with associated window 

Example 3.9. We consider the square of the Fibonacci substitution studied in Example 13.51 
q : a — > aba, b — > ab. Its inflation factor is A = r 2 , where r = 1+ 2 V ^ ■ ■ We define the valuation 
A with respect to the left eigenvector (l,r' — 1) = (1,t /_1 ) associated with the eigenvalue 
A'(= 1/r 2 ). 
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The natural decomposition (!H a ,£ft&) of its Rauzy fractal 91 is given by the solution of the 
equation 

m a = X'm a + A([a(a)] )U(X'd\ a + A([a(a)} 2 ))U(X'm b + A([a(b)}o)) 
K b = (A'$n a + A([a(a)] 1 ))U(A / ^ + A([a(6)] 1 )). 

So the previous equation can be written as 

m a = X'm a U(X'm a + T')uX'V\ b = A'9^ U (A'*K - 7- 1 ) U A'SJl 6 
m b = (A'JK« + 1) U (A'SH 6 + 1). 

The intervals 9\ a = [— l,r _1 ] and d\ b = [t _1 ,t] satisfy this equation. 

Moreover, the natural decomposition of 91 yields a tile-substitution in R, see Definition 13. 1[ 
in the following way: from (jlip we get 



(A')" 1 ^ = |J + A'-^C^OOlfc-i)) , for i G A 

(i,fc)eF l 

Note that we assume that det(M<j) = 1, thus AA' = 1, which yields 

(12) x^=(j(m j + s ji ), 

where 

(13) Sji := {AA([ffO')]fc-i) I (i> fc ) e 

Example 3.10. We continue with the square of the Fibonacci substitution studied in Exam- 
ple [33] and [3T9J Q ■ a — > aba, b — > ab. The intervals 9\ a = [—1, r _1 ] and 9^ = [r _1 , t] generate 
a dual tiling in il = R with the digit set 

_ / {0, -r} {r 2 } 
C ! {0} {r 2 } 



The tiling obtained by (fT2|) is the dual tiling of the tiling described below, where duality 
is in the sense of the cut and project scheme. Indeed, we consider as in Remark 13.41 the 
tile-substitution s associated with a with prototiles the intervals Tj of length Zj, for i G A, by 
recalling that the vector v\ = (l a , l b ) is the positive left A-eigenvector of M a such that l a = 1. 
Observe that we can associate with the two-sided fixed point u = . . . U-\.uqU\U2 ... of the 
substitution a the tiling 

According to Section [3.2[ we define the * map as: 

★ : Of A) — > R, A H> A'. 



Since 5(vw) = 5(v) + 5(w) and S(a(w)) = X5(w), for any word v,w G A*, it follows that 
A(u>) = (5(w))*, for all w G A. Hence, by ([13]), one gets 

(^)// A '(£,,),,. 

Example 3.11. One checks that {T^*j)Jj = A _1 (<Sjj)jj for the matrices of Example 13.71 and 
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Figure 3. The path associated with aaba (left), an example of a finite strand 
with coding word baab (right). 

4. Dual maps of substitutions 

In this section, we present a notion of substitution whose production rule is a formal trans- 
lation of (fT2j) . 

4.1. Generalised substitutions. We follow here the formalism introduced in [2[[3U] defined 
originally on a d-letter alphabet. We restrict ourselves here to the case d = 2. Let A be the 
finite alphabet {a, b}. 

Finite strand Let (e a ,efe) stand for the canonical basis of M 2 . One associates with each 
finite word w = w\W2 ■ ■ ■ w n on the two-letter alphabet A a path in the two-dimensional 
space, starting from and ending in A{w), with vertices in {A(w\ . . . Wi) \ i = 1 . . . n}: we 
start from 0, advance by ej if the first letter is i, and so on. For an illustration, see Figure [3] 
(left). 

More generally, we define the notion of strand by following the formalism of [6j . A finite strand 
is a subset of M 2 defined as the image by a piecewise isometric map 7: — > K 2 , where 
i,j G Z, which satisfies the following: for any integer k G [i,j), there is a letter x G {a, b} 
such that 7(fc + 1) — ~f(k) = e x . If we replace [i,j] by Z, we get the notion of bi-infinite strand. 
A strand is thus a connected union of unit segments with integer vertices which projects 
orthogonally in a one-to-one way onto the line x = y (see Figure [3]). In particular, the path 
associated with a finite word w such as defined in the previous paragraph is a finite strand. 

We introduce the following notation for elementary strands: for W G Z 2 and i G A, we set 
(W,i) = {W + Xe { I < A < 1}. 

Any bi-infinite strand defines a bi-infinite word w = (wk)kez G {a, 6} z that satisfies j(k + 
1) — 7(fe) = e Wk . The corresponding map which sends bi-infinite strands on bi-infinite words 
is called strand coding. 

This allows us to define a map on strands, coming from the word substitution, by mapping 
the strand for w to the strand for a(w). In fact, this map can be made into a linear map, in 
the following way. Let a be a substitution on A. Let us recall the notation for i G {a, 6}: 

a(i) = a{i)x ■■■er(i)| ff ( i )| = a(i)[k - l]a(i) k a(i) k+ i ■ ■ ■ a(i)\ a (i)\ ■ 

Definition 4.1. [21 [39] We let Q denote the real vector space generated by elementary strands. 
Let Ei(o~) be the linear map defined on Q by: 

E 1 {a){W,i) = Y,(M a .W + A(o-(i)[k-l}),o-(i) k ). 

k=l 
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Figure 4. The segment (0, a*) (left) and the segment (0,6*) (right). 
We call E\(a) the one-dimensional extension of a. 

Definition 4.2. Let a be a primitive substitution. The strand space X CT is the set of bi- 

infinite strands r) such that each finite substrand £ of r\ is a substrand of some E\(a) n (W, x), 
for W G I?, n £ N and x £ {a,b}. 

4.2. Dual maps. From now on, we suppose that a is a unimodular substitution. In the 
sequel we will assume that a has determinant +1. In view of Proposition 12.21 this is no 
restriction: if a has determinant —1, we will consider a 2 instead. 

We want to study the dual map E\(a) of Ei(a), as a linear map on Q. We thus denote by Q* 
the dual space of Q, i.e., the space of dual maps with finite support (that is, dual maps that 
give value to all but a finite number of the vectors of the canonical basis) . 

The space Q* has a natural basis (W, i*), for i = a,b, defined as the map that gives value 1 to 
(W,i) and to all other elements of Q. It is possible to give a geometric meaning to this dual 
space: for i = a,b, we represent the element (W,i*) as the lower unit segment perpendicular 
to the direction a of the unit square with lowest vertex W. By a slight abuse of notation, 
(W,i*) will stand both for the corresponding dual map and for the segment, i.e., 

(W, a*) = {W + Xe b | < A < 1}, (W, b*) = {W + Xe a | < A < 1}. 

For an illustration, see Figure [H Such a segment is called an elementary dual strand. 

The map E\{a) has a dual map, which is easily computed: 

Theorem 4.3. [2] Let a be a unimodular substitution. The dual map E*(a) is defined on Q* 
by 

El{a){W,i*) = {M-\W + A(a(j) k+1 - ■ ■ a(j) Hj)l ),f). 

j,k: a(j)[k]=i 

Furthermore, if r is also a unimodular substitution, then 

EI(o-ot) = EI(t)o E{((t). 
Example 4.4. We consider the square of the Fibonacci substitution. Let q: a i— > aba, b i— > ab. 
One has M e = ( ^ j J and M" 1 = ( \ ~2 ) ' We thuS get 

f EJ(e)(0,o*) = (0,a*) + (e b ,a*) + (-e a + 2e b ,b*) 
\ Et(g)(0,b*) = (e a - e b ,a*) + (0,b*). 

To be more precise, the definition of the map E\ (a) in [2j involves prefixes instead of suffixes, 
whereas, for i = a,b, the element (W, i*) is represented by the upper face perpendicular to 
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Figure 5. An example of a finite dual strand coded by the word babaa. The 
arrow at the end of the path indicates that we read the letters from left to 
right. 

the direction e*j of the unit square with lowest vertex W. Nevertheless, an easy computation 
shows that both formulas coincide. 

Dual strands. We can also define a notion of strand associated with this dual formalism. 
A finite dual strand is a subset of 1R 2 defined as the image by a piecewise isometric map 
7: — >• M 2 , where i,j £ Z, which satisfies the following: for any integer k £ there is 

a letter x £ {a, b} such that 

7(/c + 1) — j(k) = e a if x = a, ^y(k + 1) — 'y(k) = — e^, otherwise. 

Segments (W, x*), for W £ Z 2 , x £ {a, b}, are in particular dual strands. If we replace [i, j] by 
Z, we get the notion of bi-infinite dual strand. A dual strand is a connected union of segments 
with integer vertices which projects orthogonally in a one-to-one way onto the line x + y = 
(see Figure [5]). 

Any bi-infinite dual strand defines a bi-infinite word {wk)k& £ { a , b} z that satisfies for all k 

j(k + 1) — -){k) = e a if Wk = b, and ^{k + 1) — j(k) = otherwise. 

Similarly, any finite dual strand s defines a finite word w. The map tp* that sends finite dual 
strands on words in „4* is called dual coding. In particular the word coding (IF, a*) is the 
letter b, and the word coding (W, b*) is a. 

4.3. Dual strands. Before being able to define the notion of dual strand space in Section 
14.41 we need to recall several facts on the behaviour of E*(a) on finite dual strands. 

One bi-infinite dual strand plays here a particular role. Recall that v\ stands for a positive 
left eigenvector of the substitution matrix M a of the primitive substitution a associated with 
the inflation factor A. Let a be the frequency of a. We define S a as the union of segments 
(W, i*), for i = a, b, that satisfy 

0< (W,v x ) < (ei,v x ). 

One checks that S a is a bi-infinite dual strand. For more details, see e.g. [TT). One key 
property is that this bi-infinite dual strand is preserved under the action of E\(a). 

Theorem 4.5. [2] Let a be a unimodular primitive two-letter substitution with frequency a. 
The map El(a) maps any elementary dual strand of S a on a finite union of elementary dual 
strands of S a . Furthermore, if (V,i*) and (W,j*) are two distinct segments included in some 
Sa> f or a ^ (0) 1); then the intersection of their images by E^(a) is either empty, or reduced 
to a point. 
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Note that if a is not invertible, the image by E*{p) of a finite dual strand might not be 
connected. However, if a is invertible, connectedness is preserved: finite dual substrands of 
S a (i.e., connected unions of segments) are preserved under E^{a). Note that a proof different 
from the following one can be found in |17| . 

Proposition 4.6. Let a be a primitive invertible two-letter substitution. The map E*(a) 
maps every finite strand onto a finite strand. 

Proof. We first check that the generators E,L,L of the monoid of invertible two- letter sub- 
stitutions (see ([3])) map every finite dual strand made of two adjacent segments to some finite 
dual strand: 

' El(L)((x,a*) + (x,b*)) = (x,a*) + (x,b*) + (x - e a + e b ,b*) 

E\{L)((x,a*) + (x + e b ,a*)) = (x,a*) + (x - e a + e b ,b*) + (x - e a + e b ,a*) + (x - 2e a + 2e b , 
* EZ(L)((x,b*) + (x + e a ,b*)) = (x,b*) + (x + e a ,b*) 

El(L)((x,b*) + (x + e a - e b ,a*)) = (x, b*) + (x + e a , b*) + (2e a - e b , a*) 

' E{(L)((x,a*) + (x,b*)) = (x,a*) + (x,b*) + (x + e a ,b*) 

E*(L)((x, a*) + (x + e b ,a*)) = (x, a*) + (x,b*) + (x - e a + e b ,a*) + (x — e a + e b , b*) 

E{(L)(ix, b*) + (x + e a ,b*)) = (x + e a , b*) + (x + 2e a ,b*) 
, El(L)((x,b*) + (x + e a - e b ,a*)) = (x + e a , b*) + (x + 2e a - e b , a*) + (2e a - e b , b*) 

' El(E)((x, a*) + (x,b*)) = (x, a*) + (x, b*) 

Ef(E)((x,a*) + (x + e b ,a*)) = (x,b*) + (x + e a ,b*) 
Ef(E)({x,b*) + {x + e a ,b*)) = (x,a*) + (x + e b ,a*) 
Ef(E)((x,b*) + (x + e a -e b ,a*)) = (x, a*) + (x - e a + e b ,b*). 

Hence these generators map finite dual strands to connected unions of unit segments with 
integer vertices. It remains to check that these unions are indeed dual strands (i.e., that they 
can be projected orthogonally in a one-to-one way to x + y = 0). By Theorem 14.51 they all 
are substrands of S a . We deduce that the generators map finite dual strands to finite dual 
strands. 

Let us prove now that invertible substitutions map finite dual strands to finite dual strands. 
Let r be a two-letter substitution that maps every finite dual strand to a finite dual one. 
Now, if a = t o L, then we deduce from E\{a) = E^(L) o E*(t), that the map El(a) also 
maps every finite dual strand to a connected union of unit segments, and hence by Theorem 
14. 5[ to a finite dual strand. The same holds true for the other generators. We thus conclude 
by induction on the length of a decomposition on the generators E, L, L. □ 

4.4. Dual strand space and dual substitution. We now can introduce the notion of dual 
strand space for an invertible substitution. 

Definition 4.7. Let a be a primitive invertible substitution over a two-letter alphabet. The 
dual strand space X* is the set of bi-infinite dual strands rj such that each finite substrand £ 
of rj is a substrand of some E\{a) n {W,x*), for W G 1? , n G N and x G {a, b}. 

According to Proposition 14.61 since the image by E^(a) of a finite strand is a finite strand, it 
can be coded as a substitution via the dual coding ip* (introduced in Section [4.2p . 

Definition 4.8. Let a be a primitive invertible word substitution over {a, b} whose substi- 
tution matrix has determinant 1. The dual substitution a* is defined on the alphabet {a, 6} 
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as 

a*(x) = ip*(El(a)(Q,x*)) for x = a,b. 

One has L* : a \— > ba, b i— )■ b, L* : ai-> a6, 6 i— )■ b and .E* : a i— > 6, 6 i— >• a. The substitution L* is 
usually denoted as i?. (we will use it in Appendix [B]) , 

Example 4.9. Let £>: a i— > aba, b i— )> a& be the square of the Fibonacci substitution. One has 
f?* : a i— )• baa, b \-t ba. 

Remark 4.10. Let us observe that the substitution matrix of a* is the transpose of the 
substitution matrix of a. Indeed, the dual strand coded by a* (0, a* ) is located on the left of the 
dual strand coded by <r*(0, b*) (we use the fact that M a has determinant +1). Furthermore, 
one checks that a* is invertible (it suffices to check it on the generators E, L, L). Furthermore, 
the inflation factor of a* is equal to the inflation factor A of a. Furthermore, a and r are 
conjugate if and only if a* and t* are conjugate. 

This allows us to define a notion of dual frequency. 

Definition 4.11. The dual frequency a* of a is defined as the frequency of a* . 

The dual frequency is closely related to the algebraic conjugate of the frequency. Indeed, one 
has the following result. 

Theorem 4.12. Let a be a primitive two-letter invertible substitution whose substitution 
matrix has determinant 1. Let a be its frequency and a' its algebraic conjugate. One has 



Proof. According to Remark 14. 10[ the substitution matrix of a* is the transpose of the substi- 
tution matrix of a. The eigenvector of frequencies [a* , 1 — a*) of a* is thus orthogonal to the 
eigenvector (a', 1 — a') of M a . Indeed, they have distinct eigenvalues which are respectively 
A and A'. This yields a* a' + (1 — a*)(l — a') = 0, which gives the desired conclusion. □ 

Remark 4.13. There exist several codings for dual strands and dual substitutions that are 
possible [161 [TTl [9]. One could indeed code strands from right to left, or exchange the roles 
played by a and b. In particular, the coding developed by [16] yields where f is the 
substitution deduced from r by reading the letters in the images of letters by r in reverse 
order. The substitution matrices of these substitutions are either equal or transpose of the 
substitution matrix of a. Furthermore, their frequency belongs to {a* , 1 — a*}. Note also that 
a similar adequate coding can be introduced if the substitution matrix M a has determinant 
-1. 



5. Relations between distinct concepts of 'dual substitution' 

In this section we will compare the various notions of duality we have introduced so far. For 
this purpose we will introduce the appropriate equivalence relations among hulls, tiling spaces 
etc. But first we will have a closer look on the inverse of a word substitution. 
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5.1. Equivalences. If a is a primitive invertible substitution on {a,b}, then it is clear that 
cr _1 cannot be a substitution on the original letters a, 6, since for the Abelianisation would fol- 
low: M a M a -i = M a (M a )~ l = id, but this is impossible for primitive M a > and > 0. 
But we have the following: 

Proposition 5.1. Let a be an invertible substitution over {a, 6} . Assume that its substitution 
matrix M a has determinant +1. Then its inverse is a substitution on the two-letter alphabet 
{a~ l ,b}. Furthermore if a is primitive, then its inverse is also primitive. 

This is a point where we need det(Mo-) = 1. For det(M a ) = —1, things are more difficult 
(cj _1 would mix a" 1 , b with a, 6 _1 ). In such a case we consider a 2 . 

Proof. That <r -1 is indeed a substitution over {a _1 ,6} with Abelianisation M a -i can easily 
be checked by an induction on the length of a decomposition of a on the set of generators 
{E,L,L} (see ([3])). It suffices to consider L _1 and L x , they are substitutions over {a _1 ,6}. 

□ 

In Section [6] we will address the question in which cases the substitution cr _1 yields the same 
hull as a (i.e., the same set of bi-infinite words), up to renaming letters. We know already 
that is always a substitution on {a _1 ,6}. This leaves two possibilities for renaming the 
letters, namely a h-> a -1 , b H > b, or a i— >• 6 _1 , b i— >■ a. Both cases will occur: see Example 
15.41 and 15.51 below where two examples of substitutions are given whose respective inverses 
yield the same hull as the original substitutions, up to renaming letters. Thus we define the 
'inverse' substitution (that we call reciprocal) as T a cr _1 r a , where r a exchanges o with a -1 , and 
later (see Definition 16. ip we will regard a as selfdual, if its reciprocal is either conjugate to a 
itself, or to EaE, with E : a — > b, b — > a. This notion coincides with the notion of duality 
introduced in [9], see also in the same flavour [24] . 

Definition 5.2. Let a be a substitution over A = {a, b} with determinant 1, and let r a G 
Aut{F2), T a : a — > a -1 , 6—7-6. If a is invertible, then we call a := r a <7 _1 T a 7 1 the reciprocal 
substitution of a. 

Remark 5.3. If detM CT = 1, then M E M W M E = Mj. Let M = (£§). Then, since 
det M a = l, one has M a -i = ( l r p 9 ) , and the result follows from M w = ( - 1 ° ) M~ 1 ( ~ x ° ) . 
In particular, the substitutions a and <r have the same inflation factor. 

Example 5.4. The square of the Fibonacci substitution (see Example l3.5p is the substitution 
g: a — > aba, 6 — )• a6. Then, g^ 1 : a -1 — )• a _1 6, 6 — > a _1 66, and 

ff(a) = T a g" 1 T^ 1 (a) = T a g^ 1 (a^ 1 ) = T a (a~ l b) = oh. 

g(b) = T a £»" 1 r a " 1 (6) = T a g~ l (b) = T a (a _1 66) = a66. 
Thus the reciprocal substitution of g is ~g : a — > ab, b — >• a66. 

Now we can see that g and ~g are conjugate, after renaming the letters according to E : a — > 
6, 6 —7- a. Indeed, the word w = a yields a conjugation between g and EgE. 

aEgE(a)a^ 1 = abaaa^ 1 = aba = g(a), 
aE'gE^a^ 1 = abaa^ 1 = ab = g(b). 
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In particular, it follows that the hulls X Q and Xg are equal (up to renaming letters). This is 
a case where the letters change their role: in X e the letter a is more frequent, whereas in Xg 
the letter b is more frequent. 

Example 5.5. Consider a : a — > abaab, b — > ababaab. Then, <r _1 : a -1 — > ba~ 1 a~ 1 ba~ 1 , b — > 
ba~ 1 a~ 1 ba~ 1 ba~ 1 . Thus a: a — > baaba, b — > baababa. Here, w = b~ 1 a~ 1 a~ 1 b~ 1 yields a 
conjugation between a and W: 

wo{a)w~ l = b~ 1 a~ 1 a~ 1 b~ 1 baababaab = abaab = <r(a) 
wa(b)w~ l = b~ 1 a~ 1 a~ 1 b~ 1 baabababaab = ababaab = a(b) 

This is an example where the letters do not change their role. 

Now we want to compare hulls with tiling spaces and strand spaces. Thus we introduce 
suitable equivalences of these spaces. 

Definition 5.6. Two hulls X a ,X g over {a, 6} are equivalent, short X a = X e , if there is a 
letter-to-letter morphism r, such that X a = t{Xq) = {t(u) \ u E X Q }. 

In other words, X a and X Q are equivalent either if X a = X g , or if X a = E(X S ), where 
E : a — »• b, b — > a. 

Definition 5.7. Two tilings T,T' are called equivalent, short: T — T' , if they are similar, 
i.e., there are c > 0,t £ R such that cT + t = T' . Such a map is called similarity. Let 
s,s' be primitive tile-substitutions. The tiling spaces X s and are called equivalent, short: 
X s = X s /, if there is a one-to-one similarity mapping X s to X s /. 

By the analogue of Prop. 12.31 for tiling space, X s = X s / whenever there are T £ X s , T' G X s / 
such that T = T', with s, s' being primitive. 

We have seen in Section lBTTl that a word substitution a yields a one-dimensional tile-substitution 
in a canonical way: the substitution matrix M a yields the inflation factor A, and the left eigen- 
vector of A (unique up to scaling) yields the tile lengths. This, together with the order of 
the tiles in a, yields the digit set matrix D. Vice versa, a one-dimensional tile-substitution 
(where the tiles are intervals) yields a unique word substitution: just replace the tiles by 
symbols. Whenever a word substitution a is linked with a tile-substitution in this manner, 
we say a = s. If one of them (thus both) are primitive, then we write X a = X s whenever 
o~ = s. 

Using the coding maps (strand coding, dual strand coding), we can define equivalence between 
strand spaces and dual strand spaces as follows. 

Definition 5.8. Let a, a' be two unimodular primitive substitutions over a two-letter alphabet. 
We define X a = X^/ i/^(X CT /) = X a , and X a = X*, if ip*(X. a /) = X a . Furthermore, one has 

x CT ^x a , i/^(x ,) = v(x ff ). 

Note that 

X* r^> "V" f~^~> V 
a — ^cr* — sl<J* • 
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5.2. Equivalence theorem. Let us summarise these considerations. 

Remark 5.9. Let a be an invertible substitution on two letters, and let s be some tile- 
substitution derived from a as described above. Then, by definition, 

X CT = X a = & s . 
Now we can state the main theorem of this section: 

Theorem 5.10. Let a be a primitive invertible substitution on two letters with determinant 
1. Then 

X* "v ^ v "v 
<j — -A.0-* — e\.~a — A-o* — ^s* ■ 

A different proof of the equivalence between a* and El (a) can also be found in [9] . 

Proof. In order to show the equivalence of s* and E*{a), note that it does not matter on 
which lattice the paths in E\(a) are defined: all lattices of full rank in M 2 are isomorphic. 
We use here the notation of Section 13.21 So, let the underlying lattice be A =((}), ( £* ))z, 
rather than Z 2 . The substitution matrix M a acts as an automorphism on A. Indeed, since 
the projection tt\ from A to the first coordinate is (l,£x)[p) = a + (3£, and (1,£\) is a left 
eigenvector of M a , multiplication by M a in A acts as multiplication by A in 7L\£\\. 

(i,e x )M cr (i) =a(m a )(2) =\(a + pe x ). 

The same holds for £\. In the next paragraph, the important idea is that M a acts thus as an 
automorphism in Z[£], and that a stepped path considered here is in one-to-one correspondence 
with a tiling of the line. 

We proceed by translating the formal sum 

E\{a){W,i*) = Y. { M a 1 (W + A(a(j) k+1 ---a(j) wm ),f) 

j,k: cr(j)[k]=i 

into the language of digit sets and tile-substitutions, and into the internal space H. Then, 
multiplication by M" 1 in the integers Z{£\], embedded in G, is just multiplication by A -1 in G. 
And, by construction of the lattice A, multiplication by M~ in the integers Z{£\], embedded 
in H, is multiplication by (A') -1 = (A -1 ) -1 = A in H (we use the fact that a has determinant 
1). Furthermore, the term M~ 1 A(u), projected to H, reads in Z[£\] as A(u) = \u\ a + 
The formal sum above translates into 

x + T* i—)- {Xx + — tij n I n, j, such that the n-th letter in a(j) is i}, 

where tij n denotes the 'prefix' of the n-th letter, this means here: the digit dij n starred, that 
is, d*j in . Here x is the projection of W . In other words, this means 

x + If Xx + {Tf - d* jin | i = 1, . . . , m, n = 1, . . . , N} = Xx + {Tf - X>£ | i = 1, . . . , m}, 

where = N(i,j) = \{j \ o~(j) n = i}. This shows that s* and E*{a) are equivalent. 

The equivalence X w = X a * comes from the fact that a and a* have equivalent substitution 
matrices in SL(2,Z). Indeed, the substitution matrix of a* is the transpose of M CT , whereas 
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the substitution matrix of a satisfies M^M-^Me = Mj according to Remark 15 .31 By Theorem 
\7B and O the claim follows. □ 

6. Selfduality 

With respect to Definition 15.21 it is now natural to ask which substitutions and substitution 
hulls are self dual. 

Definition 6.1. Let a be a primitive two-letter substitution over {a, b} with determinant 1. 
If X a = X^, or X a = E(X^), then a is said to have a selfdual hull. 

If a ~ a or a ~ EaE, then the substitution a is said to be a selfdual substitution. 

By Theorem 15.101 this definition translates immediately to star-duals s*, and to dual maps 
of substitution a* . 

According to Remark 15.31 and Corollary 12.81 we deduce the following: 

Proposition 6.2. A substitution is selfdual if and only if its hull X a is selfdual. 

Proposition 6.3. If a is selfdual (with determinant 1), then (M CT ) _1 = Q~ 1 M CT Q for either 
Q = M Ta = ( - 1 °) , or Q = M Ta M E = (° and Mj = PM a P T , with either P = M E = 
or P = id. 

Proof. If a is selfdual, then M a = M w , or M a = MeM-^Me- Furthermore, by Remark 15 .3^ 
MeM-vMe = Mj. Hence one has either Mj = M a , or Mj = EM a E. 

Provided that detM CT = +1, one has M" 1 = ( _ s r with M a = [$.%). 

If M a = Mj, then one gets (M Ta M E )- l M a (M Ta M E ) = (l q ~ r ) = M a -i, since r = q. 

If M a = MeMJMe = ( ? s ) , then M Ta M a M Ta = ( _ p r "« ) , since p = s. □ 

The following theorem states that this necessary condition is already sufficient. 

Theorem 6.4. Let P,Q be as above (two possibilities each). If a is a two-letter primitive 
invertible substitution with det M a = 1, then the following are equivalent: 

(1) a is selfdual; 

(2) M a is of the form 

/ m k \ , ,/ / m k \ 

M m , k ={^ lm ) or M mtk =( k ^ ±1 ), 

k m 

where k > 1 divides m 2 — 1, respectively rn > 1 divides k 2 + 1; 

(3) Q- x M a Q = (Mo-) -1 ; 

(4) P T M a P=(M a f. 

Proof. It is easily seen that (2) 44> (3) 44> (4) by simple computation. Indeed, compute the 
inverse matrices: M~\ = ( m _I!L fc ) (resp. (ML r.) -1 = ( "is k )). These are obtained as 
Q^M^kQ. Obviously Q T = Q~ x . 

We know already that (1) => (3). It remains to prove that (3) => (1). Let a be invertible with 
substitution matrix of the form (2) in Theorem 16.41 Then, by Proposition 16.21 M w = M a , or 
M ? = PM a P. Then, by Theorem 12.61 a is conjugate to a, or to E(a), proving the claim. □ 
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Remark 6.5. By Theorem 15.101 the above result immediately transfers to selfdual tile- 
substitutions (on the line with two tiles) and to dual maps of substitutions. 

Remark 6.6. For every matrix M m ^ there exists a selfdual substitution having this matrix 
as substitution matrix, since the map a i— )■ M a is onto from the set of invertible two-letter 
substitutions onto GL(2,Z). 

The following result gives a necessary condition not only for a being selfdual, but — slightly 
stronger — for a ~ a*, in terms of the continued fraction expansion of the frequency a. Note 
that one can state analogous conditions for a = 1 — a* . 

Theorem 6.7. Let a be a primitive invertible substitution on two letters with frequency a. 
Let a 1 stand for the algebraic conjugate of a. The following conditions are equivalent: 

(1) a = a* 

(2) 2aa' = a + a' - 1 

(3) a = [0; 1 + m, U2, • • • , ^fcTW] or a = [0; 1, n\, ■ ■ ■ , rik], with the word n\ ■ ■ ■ rif. being a 
palindrome, i.e., (ni, • • • , n^) = (n^, • • • , ni). 

Proof. Assertion (2) comes from Theorem 14. 121 Assertion (3) comes (2) together with the 
known relation between the continued fraction expansion of a and its conjugate a 1 for a being 
quadratic (see [7] Theorem IB. 1[ (see Appendix |B1 below) . □ 



Appendix A. Sturmian words 



Sturmian words are infinite words over a binary alphabet that have exactly n + 1 subwords 
of length n for every positive integer n. Sturmian words can also be defined in a constructive 
way as follows. Let a £ (0, 1). Let T 1 = M/Z denote the one-dimensional torus. The rotation 
of angle a of T 1 is defined by R a : T 1 — > T , x i— >• x + a. For a given real number a, we 
introduce the following two partitions of T 1 : 

4 = [0,1 -a), / b = [l-a,l); 7 a = (0,l-a], J 6 = (l-a,l]. 

Tracing the one-sided (resp. two-sided) orbit of R™(q), we define two infinite (resp. bi- infinite) 
words for g G T 1 : 

a if R%(g) £ I u 



2a, e W ■ | b i fjR n(^) G / 2) 

a i£R£(g) eh, 
b if R n a (g) el 2 . 

It is well known ( [32] ) that an infinite word is a Sturmian word if and only if it is equal 
either to s Qj£) or to s a . for some irrational number a. 

The notation c a stands in all that follows for s QiQ , = This particular Sturmian word is 

called characteristic word. 

According to Theorem ll.lt invertible word substitutions on two letters are also called Sturmian 
substitutions. We recall that the set of Sturmian substitutions is a monoid, with one set of 
generators being {E,L,L} [451 129j . Moreover, we will use in Appendix [B] the following facts 
(see @D] and [29J): 
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Proposition A.l. Let a be a primitive Sturmian substitution with frequency a. There exists 
a Sturmian substitution r that is conjugate to a and satisfies r(c a ) = c a . 

If a is a Sturmian substitution that fixes some characteristic word, then a can be decomposed 
over the monoid generated by {E, L}. Conversely, if a is a primitive word substitution that 
can be decomposed over the monoid generated by {E, L}, then a is Sturmian and satisfies 
o~(ca) = Ca, with a being its frequency. 

Note that frequently Sturmian words are defined as one-sided infinite words. For our purposes 
it is rather natural to consider bi-infinite Sturmian words. 

Example A. 2. The Fibonacci sequences, that is, the sequences of the hull X a of the Fibonacci 
substitution a (see ([I])) are Sturmian words with frequency parameter a = v/ ^~ 1 . 

A detailed description of Sturmian words can be found in Chapter 2 of [29], see also Chapter 
6 in [35]. 



Appendix B. Arithmetic duality 

Let us now express the notion of duality in terms of continued fraction expansion. 

Theorem B.l. Let a be a primitive invertible substitution over {a, b} whose substitution 
matrix has determinant 1. The continued fraction expansion of the dual frequency a* of a 
satisfies the following: 

(1) if a < 1/2, then a = [0; 1 + n%, ri2, • • • , «fe, Wfc+i + with Hk+i > and n\ > 1 

• if n k+i > 1; then a* = [0; 1, n^+i, n^, • • • , n-2, n[ + n^ + \\ with k is even 

• otherwise, a* = [0; 1 + n^, n^-i, • • • , n2,n\,rij t ] with k odd 

(2) if a > 1/2, then a = [0; 1, n,2, n^, ■ ■ ■ , n^_i, + ri2\, with > 0, a* = [0; 1 + 
rife, rik-i, ■ ■ ■ , n^,ri2 + n^l, and k is even if ^ 0, k is odd otherwise. 

The description of the continued fraction expansion of a given in Theorem IB. II is due to [15]; 
such an a is a called a Sturm number; for more details, see pQ. See also in the same flavour 

M- 

Proof. We follow here the proof of Theorem 3.7 of [7] (see also the proof of Theorem 2.3.25 
of [2S]). 

Let c a (resp. c* ) stand for the characteristic word of frequency a (resp. a*), such as defined 
in Appendix[A] Without loss of generality (by possibly taking a conjugate of a by Proposition 
lA.ip . we can assume 

(14) a{c a ) = c a . 

Hence, according to Proposition IA.H a can be decomposed over the monoid {E, L}, i.e., 

a = L ni EL n2 ■■■ EL nk +\ 

with k > 1, ni,rik+i > 0, ri2, ■ ■ ■ , rife > 1. Note that the parity of k yields the sign of the 
determinant of M a . 
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Furthermore, one has L* = R and E* = E (see Example I4.9p . Hence the second assertion of 
Theorem 14.31 yields 

a* = R nk+1 E- ■ ■ ER ni . 

Moreover 

a* = EL nk+1 ■ ■ ■ EL ni E 
by using Ro E = Eo L. One deduces from Proposition IA, II that cr*(c a *) = c a *. 
Let us translate (|14p on the continued fraction expansion of a. For m > 1, let 

m '.= L m ^EL. 

One checks that for every f3 € (0, 1) 

m {cp) = c 1/{m+a) and G(cp) = c_i ■ 

We thus obtain that if 

(15) c r e mi ofl m2 o...o # mfc o L^+H^), 

with k,mi,--- ,m fc > 1, < /3 < 1, < 7 < 1, 7 = [0;^i,^2, • • • >A, ,,, ]> tnen 
P = [0;mi,m 2 , • • • ,m k ,m k+1 -\+t\,l^--- ,£i,-"]- 

With the previous notation, 

a = ni+ iofl n2 o...ofl nfc oL nfc + 1 - 1 . 
We distinguish several cases according to the values of n\ and n k+ i. 

• We first assume n\ > 0. Hence by (fl"5|) . one has a = [0; 1 + m, • • • , n k ,n k+ i + raij 
if > 0. If n k+ \ = 0, we use the fact that 

EoaoE = EL ni E ■ ■ ■ L nk = 9\6 ni ■ ■ ■ nk _ 1 L nh ~ 1 . 

Since a(c a ) = c a , then E o a o .E(ci_ a ) = ci_ a . We deduce from (fTS]) that 

1 - a = [0; l,ni, • • • ,n fc _i,n fc ]. 

Since 

1 1 1 

a 



1 _i_ l-Q 1 j \=°k 1 j 1 ' 

q 1 + l-(l-a) 1 + 3 i 

we get 

a = [0; 1 + ni,n 2 , • • • ,n k ,ni\. 

Furthermore, 

E oa*oE = L nk ^E- ■ ■ EL ni = 9 l+nk+1 6 nk ■ ■ ■ 9 n2 L n ^\ 
Since a*{c a *) = c a *, then E o a* o £7(ci_ a *) = ci_q,*. We deduce from (fT5|) that 
1 - q* = [0; 1 + n fc+1 ,n fe ,--- ,m +ra fc+1 J. 

From 

1 

a 



1-a* 1 

we obtain 

a* = [0; l,n k+1 ,n k , ■■■ ,n 2 ,ni + if n fc+i > 
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and 

a* = [0; 1 + n k , n fc _i, • • • ,n 2 ,n 1 ,n k ] if n k+1 = 0. 
• We now assume n\ = 0. One has 

a* = iT^E • • • R n2 E = SL"^ 1 • • • EX" 2 . 

— We assume n k+ \ > 0. One has a = [0; 1, n 2 , • • • , n&, n^+i]. Moreover 

a* = 9i0 nk+1 ■ ■ ■ 6 n3 L ri2 ~ l . 

We thus get 

a* = [0; l,n k+ i, ■ • -n 3 ,n 2 ]. 

— We now assume n k+ \ = 0. One has 

E oaoE = L n ' 2 E ...L nk = 6 n2+1 ■ ■ ■ nk _ 1 L nk ~ 1 . 
Since a(c a ) = c a , then E o a o E{c\- a ) = ci_ a . We deduce from (fT5|) that 
1 - a = [0;n 2 + l,n 3 , • • • ,n k _ x ,n k + n 2 \- 

Since 

1 

a = ; — , 

1+ t^T 

we get 

a = [0; l,n 2 ,n 3 , • • • ,n fc _i,n fc + n 2 J. 

Furthermore 

a* = ER nk E ■ ■ ■ R n ' 2 E = L nk E ■ ■ ■ EU 12 = 1+nfc • • • d^U 12 ' 1 . 
We thus get 

a* = [0; 1 + n k ,n k -i,--- ,n 3 ,n 2 + n k \. 
In conclusion, we have proved that 

(1) a* = [0; l,n k+1 ,n k , ■■■ ,n 2 ,m + n k+l \ with k even if a = [0; 1+nj, n 2 , • • • , n kl n k+1 + nij 
with ni > and n^+i > 0; 

(2) a* = [0; 1 + n k ,n k -i, - • • ,n 2 ,ni,n k ] with A; odd if a = [0; 1 + ni, n 2 , • • • , n^TW] with 
n\ > 0; 

(3) a* = [0; l,n k+ i, ■ ■ ■ n^, n 2 ] with k odd if a = [0; 1, ra 2 , • • • , n^, n^+i] with n^+i > 0; 

(4) a* = [0; 1+njfe, n fc _i, • • • , n 3 , n 2 + n fc J with A; even if a = [0; l,n 2 ,n 3 , • • • ,n fc _i,n fc + n 2 J. 

□ 



We now give the proof of Theorem 14.121 we want to prove that 



a 



2a' -V 

where a' is the algebraic conjugate of a. 



Proof, (of Theorem \4-12 ) Our proof is inspired by the proof of Theorem 2.3.26 in |29j . Let 



a be the frequency of a primitive two-letter substitution a. Without loss of generality, we 
assume again a(c Q ) = c a . We will use the fact that if 7 = [a±; a 2 , • • • , a n , a?], then 

-1/7' = [a n ;a n -i,- ■ ■ ,ai,o n ]. 
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• We first assume a < 1/2. According to Theorem lB.H a = [0; 1+ni, n 2 , • • • , n k , n k+ i + ni\. 
with rifc+i > and n\ > 1. Let 

7 = [n 2 ;n 3 ,-- - ,n k ,n k +i+n 1 ,n 2 \. 

We have 1/a = 1 + ri\ + I/7 and 1/a' = 1 + n\ + 1/y. One has 

-I/7' = [ nfc+1 + ni;njfc, • • • ,n 2 ,n k+1 + mj. 

We deduce that -(I/7' + "-1 + n k+ i) = [0;n k ,--- ,n 2 ,ni 

— If n^ +1 > 1, then a* = [0;l,n k+ i,n k , ■ ■ ■ , "2> ni + Wfc+ljj by Theorem lB.il We 
obtain 

CK* 1 

n fc+ i - (1/y + m + n fc+ i) = 1 - 1/a', 



1 - a* 4r - 1 

a* 

and thus 



* a'-l 

a 



2a' - 1 ' 

- If n fe+ i = 0, then a* = [0; 1 + n k ,n k -i, ■ ■ ■ ,n 2 ,ni,n k }. One gets -(1/7' + ni] 
[0; n k , ■ ■ ■ , n 2 , n\\. We deduce that 

' l/ 7 ' -m = 1 - 1/a 



l/a* - 1 
and thus 



* a'-l 

a 



2a' - 1 ' 

If a > 1/2, then a = [0; 1, n 2 , n^, • • • , n k -\,n k + n 2 ], with n k > 0. Let 

7 = N;«4, • • • ,n k _i,n k + n 2 ,n 3 \. 

We have y^ = n 2 + I/7 and = n 2 + 1/7' . One has 

-l/ 7 ' = [ nfc + n 2 ;n fc _i,--- ,n 3 ,n fc + n 2 J. 

We deduce that —(1/7' + + n 2 ) = [0;nfc_i,--- ,n3,n 2 + nfcj. 
One has a* = [0; 1 + n k ,n k _i,--- ,ri3,n 2 + n k \. Hence 



1/a* = l + n k - (1/y + n k + n 2 ) = 1 - - 

and thus 



a' 



* q'-i 

a = — ; . 

2a' - 1 



□ 
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